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INTRODUCTION 



As a result of numerous recent investigations (refer- 
ences 1-6), we are now sufficiently well acquainted, from 
an experimental viewpoint, with the laws of turbulent flow 
in smooth pir>es, channels, and along smooth plates. Laws 
of general applicability have "been found for the velocity 
distribution and the frictional resistance, and with the 
aid of these we are now in a position to handle the prob- 
lem of turbulent flow for any cases of velocities and di- 
mensions that may arise in practice. 

As far as turbulent flow in rough pipes and channels 
and along rough plates is concerned, and this type of flow 
is of greater practical interest than that along smooth 
walls, we are likewise in possession of extensive test re- 
sults accumulated within the past few years - which re- 
sults only very recently, however, have been evaluated so 
as to yield a few laws of general applicability. 

On account of the considerable importance of the prob- 
lem of surface roughness in the solution of many problems 
of engineering, the gathering of further test data on the 
frictional drag is urgently necessary. In what follows 
will be presented a new experimental method which, with the 
aid of the generally applicable laws that have already been 
found, enables surface roughness tests to be carried out in 
a very simple manner and the results to be directly applied 
to practical cases* 

In the formulation of the universal law of velocity 
distribution at a smooth or rough wall, it was found very 
useful to introduce a non-dimensional velocity cp = u/v + 
and a nondimensional distance from the wall T] = yv^/u 
* 

"Experiment elle Unt er suchungen zum Rauhigkei t sproblem. ,f 
Ingeniour-Archiv , vol. VII, no. 1, February 1936, pp. 1-34. 
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(reference 2) where u denotes the velocity, y the dis- 
tance from the wall, and v* = JTj p is a velocity comput- 
ed from the shearing stress T at the wall and may "be de- 
noted as the "shearing stress velocity, 11 p is the densi- 
ty, and v the kinematic viscosity. 

The extensive test results of J. Nikuradse (reference 
6) on the velocity distribution in smooth pipes have shown 
that the nondimens ional velocity thus formed is always the 
same function of the nondimen sional distance from the wall 
CP = cp (T| ) j the plot of cp against log T) being a 
straight line: 



o r 



cp 
u 

V * 



A + B log T| = 5.5 + 5.75 log T) 



+5.75 log — 



(smooth) (1) 



Theoretical considerations (reference 3) show that 
this straight-line law may "be expected to hold if the ef- 
fect of the viscosity on the turbulence is vanishingly 
small. It may therefore be assumed that equation (l) re- 
mains valid up to the highest attainable Reynolds Numbers. 
The velocity-distribution law (l) is a so-called boundary- 
wall law; that is, the velocity in the neighborhood of the 
wall depends, besides on the viscosity and density, only 
on the distance from the wall and on the shear stress at 
the wall, but it is not influenced by any process occur- 
ring at some great distance from the.position under consid- 
eration, in the pipe or channel; for example, not on the 
relations existing at the other wall. (This fact is taken 
into account in applying the results of pipe and channel 



Conversely, given the known values 



tests to towed plates.) 
of u, y, p, and v , the shear stress T 1 at the wall may 
be computed from the velocity-distribution law, a fact 
which we shall make use of in the carrying out of the 
tests. 



The frictional resistance law is closely connected, 
as we know, with the velocity-distribution law and is found 
to be 



1 



log (Re</M - 0.8 (smooth) 



(2) 



where X = -~ 0 is the friction factor and Re = 

dx pu 2 v 
is the Reynolds Number, This law has the same range of 
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application as the velocity-distribution law (l) and there- 
fore likewise applies to arbitrarily high Reynolds Numbers. 

If we plot similarly the relation cp = f (log T) ) for 
rough pipes, then from the measurements of Nikuradse (ref- 
erence 10) on pipes artificially roughened with sand, a 
straight line it obtained for the velocity distribution 
for each value of relative roughness k/r (k = absolute 
roughness, r = pipe radius) and each Reynolds Number. It 
is found, in particular, that all the straight lines run 
parallel, the common slope B = 5,75 being the same as for 
the smooth pipe. Using the identity 

yv* y v*k 
V k u 

the universal velocity-distribution law for rough pipes 
may therefore also be written in the form 

m * A + S.fS log Z, A = A (~) (rough) (3) 

^ k v V / 

According to the measurements of Nikuradse the magnitude 
of A depends only on the number V*k/tf , involving the 

"shear stress velocity" v* , the roughness k, and the 
kinematic viscosity V. From equation (8) the resistance 
law for a rough wall is obtained a s 



A = (2 log Z + a)* -2 (rough) (4) 

where a similarly depends on v*k/u. 

The velocity-distribution and resistance laws (3) and 
(4) for rough walls were extensively investigated by 
Nikuradse, whose object it was to investigate a certain 
simple type of roughness over a very large range of Reynolds 
Numbers. The roughness chosen was sand of grain sizes^ 
fc s = 0.8* 0.4, 6.2, and 0.1 millimeter, the sand being 
glued on as densely as possible to the inner surface of 
brass pipes of various diameters (d m 2.5, 5, and 10 cm) by 
means of" lacquer. The greatest relative roughness was 
therefore, r/k s = 15, and the smallest, r/k s = 500. 
Figure 1 shows the dependence thus found of the resistance 

u d 

or friction factor A on the Reynolds Number Re = — 

and the relative roughness r/k g , X = f (Ho, r/k g ), for a 
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range of HkfnoX&s Numbers from 10 3 to 10 6 . It may "be seen 
that for the friction factor three different regions may be 
made out which may "be characterized as follows: 

In region I, that for small Reynolds Numbers, 
the resistance of the rough pipe is equal to that of 
the smooth one. The roughness elements lie entirely 
within the very thin laminar layer which in turbulent 
flow adheres quite close to the wall (laminar sub- 
layer } . 

In region II (region of transition) the resistance 
increases with increasing Reynolds Number. The thick- 
ness of the laminar layer is of the same order of mag- 
nitude as the height of the roughness, so that with 
increasing Reynolds Number (that is, with decreasing 
thickness of the laminar layer) more of the roughness 
excrescences protrude through the layer and therefore 
the resistance strongly increases. 

In region III the resistance factor is independ- 
ent of the Reynolds Number and depends only on the 
relative roughness. All the roughness elements pro- 
trude through the laminar laygr and the flow may now 
be considered entirely as turbulent flow. For A we 
then have the simple formula 

A - fl*74 + 2 log ~- N ! " (5 ) 

K s J 

The velocity-distribution law (^) likewise assumes a 
particularly simple form for region III since A also 
becomes independent of v^k/u. For the special case 
of the sand roughness of ITikuradse, A - A s = 8.48 
in this region. 

The three regions thus made out apply not only to the 
sand roughness but also to many other kinds of roughness 
met with in practice. The region III where the flow may 
be considered as entirely turbulent is the most important, 
and our investigation will' be mainly concerned with this 
r egi on. 

The universally applicable laws for smooth and rough 
walls we have just considered, enable us to carry out 
roughness tests in a particularly simple manner as will be 
described below. The idea suggests itself of establishing 
a "standard roughness" to which any arbitrary roughness may 
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fee referred for comparison. Such a standard would natur- 
ally "be quite arbitrary and the choice would depend on 
practical considerations. We shall not, however, consider 
this subject further here. Without setting up any stand- 
ard roughness the problem in any practical case where rough- 
ness must he taken into account may "be considered as com- 
pletely solved if it is possible to indicate the resist- 
ance for any kind of roughness "both for the case of pipe 
and channel flow and for a plate towed in a region of un- 
confined flow - the latter case "being of fundamental impor- 
tance in the study of ship hull design. In "both cases a 
necessary requirement is that it should "be possible to ex- 
tend the results directly to pipes and channels or plates 
of other dimensions. It should he possible for any given 
geometric form of the roughness, to compute the frictional 
resistance theoretically, although such computation would 
rarely succeed in giving the required accuracy. It is all 
the more important therefore that we he in possession of a 
simple experimental method whereby, without any undue time 
expenditure, any practical problem involving surface rough- 
ness that may occur in practice, may be investigated to or- 
der, so to speak. An apparatus of a simple type has been 
developed in Gottingen and has the form of a rectangular 
channel with three smooth walls and one interchangeable 
rough wall. This apparatus is found to be particularly 
useful for ship-design problems, since by a simple exten- 
sion of the results obtained in this channel, they may be 
directly applied to the case of a towed plate. In this 
manner the* plat e-t owing tests that in many cases involve 
elaborate testing and much time and money expenditure, may 
Be dispensed with • 

In order to have the characteristics of the new exper- 
imental procedure clearly brought out and at the same time 
to obtain some further results on the nature of florin 
rough pipes and channels, we have investigated a series of 
uniform, geometrically similar, simple types of roughness. 

In the tests of Hikuradse the roughness was charac- 
terized by a single parameter, namely, the absolute rough- 
ness k = k s = size of sand grain. On account of the va- 
ried character of the types of roughness that may come up 
in practice, it will not be found possible, in most cases, 
to use a single parameter but at least one more will be 
found necessary, namely, the roughness density, which is 
the number of individual roughness elements per unit of 
area. (In the Hikuradse tests the roughness density had a 
constant value near' the maximum.) It is therefore also 
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necessary to investigate how, for a similar, geomet r i cal 
form, the frictional resistance depends on the density of 
the roughness., and this question we shall likewise consider 
in the tests described "below. 



I. DESCRIPTION 0? TESTS 
Test Set-Up 



The tests were conducted on the test set-up shown in 
figures 2, !5, and 4, This is essentially the same as that 
previously employed "by Nikuradse in his tests on smooth 
and rough pipes. A detailed description of the entire set- 
up is given in the V . D . I • , For s chung sh eft 256 (reference 6), 
so that we need here only consider wherein our apparatus 
differs from his. 

The centrifugal pump fcp delivers the water from the 
stored water channel vk into the water tank wk. From 
the latter the water flows through the pipe zr and the 
starting run channel ak , and reaches the channel mk 
where the measurements are taken and then flows through the 
velocity-measuring apparatus gm hack into the stored 
water channel vk. The channel on which measurements are 
taken is more clearly shown in figure 4. It consists of 
two parts, each 3* 20 meters long, screwed to each other at 
the flanges. Both channels have the same cross section 
4 # 0 by 17.0 centimeters. The actual measurements are made 
at the rear portion, the forward portion serving as a take- 
off run. The walls of the latter portion are smooth. This 
take-eff channel, which originally was intended as a meas- 
uring channel, was welded together of U "beams and sheet 
iron. The latter method of construction turned out to he 
too inaccurate, however, since, for reliable results to he 
obtained, it is necessary that the height of the channel 
along its entire length should net vary by more than l/lO 
millimeter (see below). The second channel was there- 
fore constructed of cast steel and consists of two parts: 
1 - a rectangular open run with smooth walls and grooves 
at the sides, in which the test plate is fitted; and 2 - 
a cover which is placed above the test plate and is tight- 
ly screwed onto the lower half of the channel (fig. 4). 
To obtain sufficient tightness between cover and test 
plate a linen packing sprinkled with red oxide of lead was 
placed between cover and test plate* By careful workman- 
ship of the walls and side grooves at the lower half of 
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the channel, uniformity in the height of the channel all 
along the channel length was achieved. As a protection 
against rust and to attain sufficient smoothness of the 
walls, the inside of the channel was spread over with 
"inertol" and rubbed smooth with sandpaper. The very 
slight roughnesses that still remained do not enter into 
the results as errors "but are taken into account in the 
evaluation of the results "by the method described "below. 
The wall thus treated may practically "be considered as 
hydrauli ce lly smooth. 

At eight positions along the smooth wall, 40 centime- 
ters apart, the statical pressure was taken. At each sec- 
tion three orifices were "bor.ed so as to "balance out any 
small errors at each section "by taking an average. At each 
section at which the static pressure was taken, the channel 
was provided with small windows through which the flow 
could "be observed, It was also possible to insert a pitot 
tube through an orifice in the smooth wall so that the ve- 
locity distribution could be measured at various sections 
of the channel. With the aid of a block and pulley ar- 
ranged above the channel cover, the latter could easily 
be lifted up after unloosening the screws, and the test 
plate interchanged with another. 

Types of Roughness Elements Used 

The test plates upon which the roughness was to be in- 
vestigated were of such dimensions as to fit accurately 
into the channel, their lengths being 320 centimeters and 
width 17*8 centimeters, of which 5 millimeters were allowed 
on each side for fitting into the side grooves. These 
plates were of sheet metal 5 millimeters thick, on which 
was placed square bar iron of 15 millimeters thickness, so 
that the total plate thickness was 20 millimeters. Care 
was taken to see that the plates were very even. 

Table I gives the data on all the rough plates inves- 
tigated, and figure 6 shows the appearance of each plate. 
Six groups of roughness elements were used, and each group 
at several roughness densities. The dimensions and arrange- 
ment of these roughness elements are indicated on figure 5; 
I, spheres, diameter d = 0.41 cm; II, spheres, diameter 
d = 0.21; III, spherical segments; IV, cones; V, "short 
angles"; VI, "long angles"; and in addition, a plate with 
"Hamburg sand" of a mean grain size k = 1.35. This is the 
same sand that was employed by G-. Kempf (reference 17) in 
his towing tests on rough plates. In the preliminary tests 
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the same kind of sand was employed as that used by ITikuradse 
in his tests on pipes. To prepare a sand plate the plate 
which was at first well smoothed, was sprayed* with "iner- 
tol" and the sand was spread directly over it, the lower 
sand layer strongly adhering to the plate. After shaking 
off the excess sand the plate was allowed to dry for one 
day and then it was again sprayed with inertol. After a 
second drying the sand layer is found to adhere strongly 
enough to the plate so that it is not torn loose "by the 
water streaming through. 

The other roughness elements were all soldered onto 
the plates. Fo'r this purpose the plates were first cov- 
ered over with a thin layer of tin. The spheres, sphere 
segments and cones, except when they were most closely 
packed together, were soldered on using sheet aluminum 
patterns, the short and long angles and the closely packed 
spheres and spherical segments "being soldered on without 
any pattern. The spheres could he obtained as round shot. 
Before "being soldered on they were given a thin tin coating. 
The spherical segments were stamped out of galvanized sheet 
iron 0,? millimeter thick and the cones were similarly 
stamped out of thick sheet "brass. The angles were formed 
out of galvanized sheet iron 0#3 millimeter thick. Ylith 
the exception of the plate having the long angles, the el- 
ements were distributed on all the plates as shown in fig- 
ure 5. The distance between any two elements of the same 
horizontal row is the same as the distance between two hor- 
izontal rows. The rows are staggered as shown in the fig- 
ure so that each element is centered "between the others. 

Pr el iminary Tests 

In the preliminary tests there was only the channel 
shown on figure 4 at the left. The pressure-drop measure- 
ments made with the orifices and also with the static tubes, 
"both sides of the wall "being smooth, showed a very non- 
uniform pressure variation along the channel length. It 
was immediately recognized that the reason for this was 
the lack of uniformity of the cross sections at the vari- 
ous points of the channel. With the relatively large cross 
section and therefore small pressure drop, even very small 
variations in the cross section 1 show greatly disturbing ef- 
fects on the pressure-drop measurements as may he seen from 
the following- approximate calculation. 
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TABLE I 

Dimensions of the roughness elements k, k f , d, c, D, D/d ; the roughness 
density F r /P, F x /f and the mean channel height b v for the various 
test plates. For the meaning of d, D, k, k ! , c see figure 5. b = 
rr.epn height of rough channel; F = plate area; F r b projection of 
roughness elements on plane normal to direction of flow; F A = area of 



smooth parts of plate Detween roughness elements. 



Type of . 
roughness ! 


Plate 
No. 

J 


d 

err. 


j 

D 

cm 


D 

d 


1_ 

K 

cm 


j short angle 

I Innp* «n<p*lp 
err? k cc/iioa 


h 

! 

— 'i 


T7 1 

F r 

F 


F 


Spheres 


xil i 




0.41 


—< \ 


9.75 


0 41 




3.99 


0.00785 


0.992 


Spheres 


XI la 1.0 


io ; 


10 


1.0 




_ 


.00785 


.992 1 


Spheres 


III 


.41 


2 j 


4.88 


.41 


** 


3.99 


.0314 


.969 


Spheres 


I ! 


.41 


1 1 


8,44 


.41 




3.96 


.126 


.874 


Spheres 


II | 


.41 


0.6! 


1 .46 


.41 




3.88 ' 


.349 


. 651 


Spheres 


V! 


►41 packed tightly 


.41 




3.68 


.907 


.093 


Spheres 


VI 


.21 


1 


4.86 


.21 




3.99 


.0314 


.969 


Spheres 


IV, 


.21 


0.5 , 


c.43 


.21 
*• 




3.97 


.126 


.874 


Spherical 






I 














segments 


XIII 


.6 


4 ! 

1 


5 


.26 


'■ 1 * 


3.99 


.0087 


.969 


Spherical 




















segments 


XIV 


.8 


3 | 


3.75 


.26 




3.99 


.0155 


.944 


Spherical 




















segments'; XV 


.8 




2.5 




_ 


3.98 


.0348 


.674 


Spherical 


1 




i 














segments 


XIX 


.8 


pa eked 


tightly 


.26 




3.85 


.251 


.093 


Cones 


XXIII 


.8 


4 ! 


5 








01 On 


9 69 


Cones 


, XXIV 


.8 


3 


3,75 


.375 


.425 


3.98 


.0189 


.944 


Cones 


1 XXV 


.8 




i.5 


.375 


.425 


3.95 


.0425 


.874 


Short 




















angles 


XVI 


_ 


4 




.30 


.8 


4.0 


.0151 


.998 


Short 




















angles 


| XVIII 




m 




.30 


.8 


4.0 


.0269 


.996 


Short 




















angles 


XVII 




2 


m 


.30 


.8 


3.99 


.0605 


.994 


Long 


















.995 


angles 


XX 






— 


.32 


17 


3.90 


.0538 


Long 
angles 


XXI 




4 




.31 


17 


3.96 


.0776 


.992 


Long 
















.152 




angles 


i XXII 




2 




.30 


17 


3.96 


.985 


Hamburg 














3.87 






sand 


IX 

1 


m 






.135 









1 Flate XI la was measured in larsre tunnel only. 
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The mean velocity of the flow is at most u = 8 me- 
ters per second, the hydraulic diameter 

g ;=s 4_x_crosirjectional,area = 6.4 8 cm 
h wetted perimeter 

the 'kinematic viscosity v = 0,012 cm 2 per second, so that 
the Reynolds Uumher Re ^ u | = 4.3.x 10 5 . With this 
Reynolds Number -the 1 resistance factor of the smooth pipe 

is x = -~ = 0,018 where q is the dynamic pressure 
dx q 

corresDonding to the mean velocity. The frictional pres- 
sure drop between t¥o points of measurement at a distance 
I = 40 cm apart is therefore 

*» ■ 1 s ■ * 4 « - °- u ^ 

If the channel height b is now assumed to change "between 
the two points "by the amount- Ah, the width remaining 
constant", the difference in pressure thus produced is ac- 
cording to the Bernoulli equation 

Ap' = p u Au = p u 2 ^ = 2 A£ q .. 

If the error in the ■oresa'ure drop is not to exceed 2 per- 
cent, then 

.. M- < 0.02 



or 



A p 

2 ^< 0.0022, Ah < 0.0011 h 

n 



D 

or finally, with h = 4 cm 



Ah < 0.0044 cm ft mm 



mm 



The channel height must therefore he accurate to 1/20 
throughout the channel length if the static pressure meas- 
urements arc to he considered reliable. It was not pos- 
sible to attain this accuracy with the channel as first 
constructed of welded parts and therefore the second chan- 
nel was made of cast steel as has previously been mentioned 
and it was then possible by carefully working over the 
walls to attain the required accuracy. 
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In the preliminary experiments it was also checked 
to see whether the flow was uniformly established after 
reaching the end of the first channel. The starting run 
is shorter at a rough wall than at a smooth one and it 
was found that for a small height of roughness the flow 
could not yet "be considered as uniform after the end of 
its run through the first channel so that after the second 
oast— steel channel was completed the first was used as an 
initial runway and this arrangement was kept in the suc- 
ceeding tests. It was also "brought out by the preliminary 
tests whether in making the pressure drop and velocity- 
distribution measurements in the second channel, it made 
any difference whether two smooth walls or one smooth and 
one rough wa.ll were used. Since this was not the case, two 
smooth walls were used with the initial run channel for all 
the main tests. 

The pressure drop in the measuring channel proper with 
"both walls smooth was likewise measured in the preliminary 
tests* A friction factor X, referred to the hydraulic 
diameter, was obtained that was found to "be in good agree-* 
ment with the law of resistance of smooth pipes. 

It was also investigated "by the preliminary tests with 
a very small channel of rectangular cross section (1.05 x 
5 cm) having one rough and one smooth wall whether the ve- 
locity profiles at the rough and smooth wall, respect ively , 
had any effect on each other. In evaluating the results "by 
the method described in a subsequent section (p. 15 ) , it 
has been assumed that the friction layers are formed at 
the smooth and rough wall independently of each other and 
are the same as for a channel with all smooth or all rough 
walls. This assumption may be tested in a simple way by 
plotting the velocity at each wall against the logarithms 
of the corresponding. distances from the walls. Since, ac- 
cording to the generally applicable velocity-distribution 
law the velocity is proportional to the logarithm of the 
distance from the wall for both rough and smooth walls 
the above plot must show a triangular velocity distribu- 
tion. Figure 7 shows how well it actually does give this 
di st ribution. 

Before starting the main tests there was also obtained 
for the outlet section of the measuring channel a diagram 
showing the lines of equal velocity (fig. 8). The width/ 
depth ratio of 4.25: 1 for the rectangle was thus found to 
be sufficiently large for the flow to be considered as two- 
dimensional at the center of the section to a sufficient 
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degree of accuracy, a necessary underlying assumption for 
the method employed of evaluating the results. 

Carrying Out of Main Tests 

In the main tests measurements were made to obtain: 

1. The velocity profile parallel to the short sides 

at the center portion of the outlet section. 

2. The pressure drop by means of orifices. 
Zm The temperature. 



4. The volume discharged. 



In measuring the velocity distribution the total pressure 
given by the pitot tube was calibrated against the static 
pressure at the outlet section so that the dynamic pres- 
sure was measured directly. In the pressure-drop measure- 
ments all possible combinations Pa P3 i P3 ** P4 » ■ 
P2 m P7 1 P3 "* P 4 j P 6 - P 7 were measured so that by 

talcing a mean value it would be possible to balance out 
small errors occurring in the individual pressure measure- 
ments. The pressure at point j), was not included since 

the flow did not assume a steady condition at this point. 
The arrangement of the roughness elements was such that 
the last row of elements was directly along the edge of 
the plate. The pitot tube moved in a plane l/2 millimeter 
behind the plate edge. The profile measurements were ob- 
tained, as far as possible, at the center of the space' 
between two longitudinal rows of elements. A velocity pro- 
file was also measured, however, at the center behind a 
roughness element. For most of the plates both profiles 
differed only in the immediate neighborhood of the wall. 
Only in the case of plate XII (spheres of diameter d = 4 
mm spaced at distance D = 40 mm apart) were any spe- 
cial phenomena observed in connection with the velocity 
profiles, and these we shall consider later. 



In order to keep the motor speed constant the tests 
were partly conducted using the automatic speed regulator 
described in the V.D.I., For schungshef t . 356 , partly without 
using it. In the latter case, for the purpose of main- 
taining uniformity with respect to time of the velocity 
through the channel between two measuring points, the pres- 
sure drop along the channel was controlled by a second 
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manometer and the throttle of the velocity-measuring ap- 
paratus somewhat adjusted. 



II. EVALUATION OF TEST RESULTS 
Some General Considerations 



It is the ohject of our tests to "be able to specify 
for each type of roughness investigated, a characteristic 
number "by whose aid it "becomes possible to predict the re- 
sistance for the same type of roughness at other Reynolds 

__ di_ k 
Numbers Re = u -~ and relative roughness ratios — 

(where r^ is the hydraulic radius) than those directly 
measured. The number k denotes any useful measure of the 
absolute roughness and for uniform roughness elements, is 
most conveniently chosen as the maximum height of the 
roughnes s. 

Since in the case of sand roughness (grain size denot- 
ed "by k g ) we possess a type for which the various resist- 
ance relations are known over a wide range of Reynolds Num- 
bers and relative roughness, it appeared advantageous to 
us to evaluate our results in such a manner that they could 
easily be expressible in terms of the results obtained by 
Nikuradse for sand roughness without, however, necessarily 
taking the latter as a standard type of roughness. One ob- 
jection to using sand roughness as a standard is that this 
type is not satisfactorily reproducible. If, for example, 
instead of employing lacquer for gluing the sand to the 
plates, as was done by Nikuradse, we used some other bind- 
ing material, and if the time taken for drying was chosen 
somewhat differently, another sand-grain density and there- 
fore a different resistance, would result. 

The above consideration is of little account as far as 
our method of evaluating the results is concerned, since 
we consider only the formula giving the resistance factor 

as a function of the Reynolds Number Re m u -Jl and the 

ratio ^s/ r h "but ^se no comparison values on the rough- 
ness densities and absolute roughness. The results of the 
resistance measurements of Nikuradse are given in figure 1 
where A is shown as a function of Ud/v and k s /r for a 
very large range of Reynolds Numbers and relative rough- 
ness rat io s : 
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The function X = f (Re, k s /r) assumes a particularly 
simple form in region I (see Introduction, p. 4) of hy- 
draulically smooth flow, and in region III, of completely 
rough flow. For region I, X ceases to depend altogether 
on k/r and the simple resistance equation (2) for which, 
in the case of smooth pipes, Nikuradse gives the approxi- 
mate formula 

X . 0.0032 * °; g ?* v (2a) 
Re 0 * 23 



becomes valid. In region III, the dependence on the rela- 
tive roughness may "be expressed "by the simple formula given 
by von Karman (reference 3) and "by Prandtl and Nikuradse 
(references 2 and 10), already indicated in the Introduc- 
tion, namely, 

X = (l.74 4- 2 log ^ 2 (5) 

the graph of which is shown in figure 9 where r^ denoted 
the hydraulic radius = 2 X cross section/wetted perimeter. 

For the region of transition XI, where X depends on 
Re as well as on r/k, we have as yet no corresponding 
analytical formula for X. 

The region of greatest practical importance is "by far 
region III, which extends over a large range of Reynolds 
Numbers and not too small values of absolute roughness. 
For this region the relations fortunately assume an excep- 
tionally simple form. According to the measurements of 
Nikuradse on sand roughness, the purely rough type of flow 
occurs for the conditions 



> 70 or Ufc > 125 2* 

In our own series of measurements the relative roughness 
k/r-^ and the Reynolds Numbers were sufficiently large so 

that the square law of resistance could "be held to apply 
throughout • 
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Explanation of the Method of Evaluating the Results 

Our problem, namely, to specify for each of the types 
of roughness elements investigated, the resistance of those 
of similar geometric form hut other relative roughness - 
that is, to extend the results to pipes and channels of 
other diameters, may he considered as solved if, for each 
type of roughness we can indicate the value of the "equiv- 
alent sand roughness' 1 k , that is, the size of sand 
grain of Hikuradse that has the same resistance as the 
roughness inve st igat ed. No particular physical signifi- 
cance is to he attached to this equivalent roughness hut 
it is rather to he considered as a convenient magnitude 
that allows us, with the aid of equation (5), to extend the 
results immediately to the computation of the value of A 
for other types of roughness. Instead of k g we might 
also use the nondimensional factor 



where k denotes the actual ahsolute roughness (approxi- 
mately the maximum height of roughness) of the type under 
consideration. 

This process of reducing any arhitrary roughness com- 
putation to that of the equivalent sand roughness may he 
made clearer hy an example. We shall take the measure- 
ments given hy Hopf (reference 7) and Eromm (reference 8) 
on three different roughnesses, namely, wire netting, "waf- 
fle" sheet metal, and saw profile. The measurements were 
taken for different channel heights, that is, for differ- 
ent relative roughness. For the entirely rough or turbu- 
lent type of flow the straight-line equation of slope 2 



must apply. This equation is found to he fulfilled in a 
satisfactory manner and the values derived for the constant 
a are given in tahlc II. By comparison of equation (5) 
with (3) there is ohtained for the reduction coefficient 



k 




(3) 



a = k ^ / k the relation: 



2 lo 



g a = X # ?4 m a 



The corn-out ed values of a and 
table II, 



are likewise given in 
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TABLE II 



Reduction of Measurements of Hopf and Eromm 
to Equivalent Sand Roughness 



k = absolute roughness, 


i 




- 2 log 

k 


k ' K a 


= equiva 


jl en t 


sand 


roughne 


s 


c 1 


a = 1 


s _ 




Type of roughness 


I 


k cm 




a 


a 


k s cm 


Wire netting 


h 


0.0115 






0.96 


2.46 


0.028 


TTaf fie- sh eet metal 




.0427 






1 . 36 


1.52 


.065 


S aw profile 




. 15 






1.48 


1. 34 


. 201 



i 



In this manner it is possible to reduce every rough- 
ness within the range of the turbulent flow to the equiva- 
lent sand roughness if the resistance coefficient is known 
as a function of the relative roughness. 

To determine the equivalent sand roughness k s it is 
not necessary, however, to mea.su.re every roughness at vari- 
ous relative roughness ratios, a single measur emen t at one 
relative roughness only "being sufficient as "becomes immedi- 
ately evident from what follows. The amount of testing 
required is thereby greatly reduced. 

In our tests with the rectangular channel having a 
smooth and rough wall, the determination of k s is some- 
what complicated by the fact that from the pressure-drop 
measurements we know only the over-all channel resistance 
which is made up additively of the rough and smooth wall 
resistances. We therefore still require further informa- 
tion on the resistance of the smooth wall in order to be 
able to compute the rough wall resistance which alone in- 
terests us. The resistance of the smooth wall is obtained, 
as will be explained directly below, with the aid of the 
generally applicable velocity-distribution law given above 
from the measured velocity profile in the vicinity of the 
smooth wall. 

Denoting by T r and Tg the shear stress at the 
rough and smooth 'wall , respectively, there is obtained 
from the equilibrium of shear stresses at the wail with 
the pressure drop, a.s a first equation for the determina-* 
tion of the unknowns T and t^. the relation 
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where "b denotes the height of the channel. The velocity 
distribution lav; for the friction layer at a smooth and ro 
rough wall, respectively, is, according to Prandtl (see 
Intro duct ion ) 

v rj.y 

= 5.5 + 5.75 log — — (smooth trail) (7) 

JL = A + 5 . 75 log y (rough wall) (8) 
v *r - 

where v end v denote the ,r shearing stress veloci- 
*g *r 

ties" at the rough and smooth wall, r espect i vely , namely, 

^ v *e - ff- r ** = rf 

and A for the type of roughness in question is a charac- 
teristic function which is constant for the region of rough 
flow hut depends on the roughness coefficient v*k/i; in 
the region of transition from smooth to rough flow. For 
sand we have for the rough flow, according to the results 
of Mkuradsc, the value of A«, = 8.48 and hence the veloc- 
i ty-di strihution law for thr case of sand roughness "becomes 

= 8.48 + 5.75 log XT ( 8a ) 
v *r k s 

Tie now make use of the velocity-distr ihut ion law for 
smooth vails (?) to determine the shear stress at the smooth 
wall from the measured velocity profile. Plotting the meas- 
ured velocities against the logarithms of the distances 
from the smooth wall, there is obtained a straight line 



u - nig + ng log y 

(7a) 



whose slope n g immediately gives Tg« Comparison of (7) 
with (7a) ^ives n~ = 5.75 v.^ 
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The slope of the profile at the smooth wall can he oh- 
tained graphically to a sufficient degree of accuracy. In 
this method of determining v m it is assumed that the 
friction layers at the smooth and rough wall have no mutual 
effect on each other hut that they have the same form as 
in the symmetrical tunnel. That this is true was shown in 
the preliminary tests (fig. 7). We may, instead of the 
shear stresses themselves, use the corresponding velocities 



v 



and v and equation (6) then takes the form 



* r * g 



or 



v a+ Y a = l *p (6a) 

* r *g p dx 



/."b dp a pb dp / n g \ 2 , ft v 



This is a first method for the determination of the unknown 
merits. 



v *r upon which was "based the evaluation of our measure- 



Before we proceed to determine from the value of 
the coefficient a = for reducing to equivalent rough- 

ness, we shall indicate yet a second and third method for 
the determination of v* r vrhich are independent of the 
first one just explained and may therefore he used as a 
check . 

In the same manner that we determined v^g with the 

aid of the v el o c i ty- di s t r i hut i on law for a smooth plate 
(7), we may use the corresponding law for rough plates (8) 
to determine v* r * Plotting the measured velocity against 
the logarithm of the distance h - y from the rough wall 
we similarly ohtain a straight line 

u = m r + n r log (h - y) 

from whose slope n r , hy comparison with (8), we find 

* r 2 5.75 

The slope n r may also he determined very simply hy the 
graphical method. We have thus determined v* r from 

equation (10) and v *r 2 fron equation (11) and "both val- 
ues must agree. 
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It was found that the value of v *r 2 generally ex- 
ceeded somewhat the value of v ^t 1 C^y a"bout 5 percent) 

and this may he explained as follows. In the measurements 
of Nikuradse on smooth and rough pipes the result was ob- 
tained that the velocity profiles of those portions lying 
more toward the center of the pipe generally showed a 
slight systematic deviation froffi the velocity distrihution 
equations (7) and (8). (3ee 7. D.I,, For schungshef t 356, 
fig. 24, and V.D.I*, Fo r s chungshef t 361, fig. 14.) Most 
of the measured profiles ran somewhat steeper than the 
straight line with slope 5.75. In determining v *g and 

v r from equations (9) and (11), respectively, we should 

have taken a somewhat larger value than 5.75, the precise 
value, however, being unknown. Our obtained values for 
and v. v are therefore consistently somewhat too 

large, amd for this reason ^^ Vl is found according t* 
(10) to be too small. This explains why v *r 2 was gener- 
ally obtained larger than v* r _^ . The actual value of v* r 
must thus lie between v^~ and v^^ and we therefore 
took the mean of these two values 

v*r = \ (t* Pi + ^*r 2 ) ( 12) 

as being'the nearest approach to the actual value. 

Finally, a third method for determining v* r is 
given by the position of the maximum of the nonsymmetrical 
velocity profile. If b x and b 2 denote the distances, 
respectively, of this maximum from the smooth and rough 
wall (b x + b 2 = b ) , then, since in completely turbulent 
flow the shearing stress is a linear function of the dis- 
tance from the wall and at the maximum velocity t - 0, 



or. 




If v^g has been determined from the velocity-distribu- 
tion law, then v* r may be computed from the above equa- 
tion. In evaluating our results we made use of this meth- 
od for the determination of v r only as an incidental 
che ck. 
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We now proceed to consider hew the results obtained 
with our types of roughness are reduced to the equiva- 
lent sand toughness of Nikuradse. For this piirpose we 
determine, for each type of roughness, the constant A of 
the universal velocity distribution law. 2 From equation (8) 

A * ~ 5.75 log I (13) 

v *r K 

where k may- denote the maximum height, of -the roughness 
under consideration. The numerical values of k are 
given in table I. In equation (13) setting k = k^/co, 
we find by comparing (13) with (8a) the expression for a: 

5.75 log a = 8.48 - A (14) 

connecting the conversion factor a with the constant A 
of the velocity-distribution law. 

In evaluating our test results we computed, for each 
velocity profile, a mean value of A from equation (13) 
using for v^ r the mean of the an ^- v *r 0 values 

derived above. With the relatively large roughness ele- 
ments employed the flow obtained was of the purely rough 
type throughout. This may be seen from the fact that the 
values of A for each of the profiles differ only slight- 
ly from those for a rough plate. (See table III, p. 40.) 
Out of these values a mean value of A was obtained from 
all the measurements on one plate (table IV) and with this 
value of A the coefficient a was determined. 

In evaluating the results of the measurements, par- 
ticular consideration must be given to the measurement of 
the distance from the rough wall. This distance shall be 
so defined as to be equal to the distance from a hypo- 
thetical smooth wall that replaces the rough wall in such 
a rnanner as to keep the fluid volume the same. With both 
walls smooth the volume of the channel is accurately known 
(cross section 4.00 by 17.00 cm). Since the volume of the 
roughness elements may be computed, the mean cross section 
and therefore the mean channel height with rough and smooth 
wall is known. These mean channel heights are given in 
table I for each of the roughness types. 



This is identical with the function X employed 
by Prandtl (reference P). 
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Test Results 

The results of cur tests are summarized in figures 10 
to 15 and table III. The measured velocity profiles were 
evaluated by the method described above. In figures 10 to 
15 the results are plotted for each of the types of rough- 
ness (spheres, spherical segments, cones, short and long 
angles) according to the universal velocity-di stribut ion 
law (u/v #r against log y/k). All profiles measured for 
the same plate must lie on a single straight line insofar 
as completely turbulent flow is assumed to hold. The fig- 
ures show that they actually do so to a sufficient degree 
of accuracy, each rough plate corresponding to one such 
straight line and all straight lines running parallel as 
required by the veloci ty-di stribut ion law. The slope of 
these lines is, on the average, somewhat larger than the 
value 5.75 of the velocity-distribution law. This fact, 
and the manner in which it has bem taken into account, has 
already "been discussed cbove. That the profiles obtained 
"by us do not coincide with the straight lines of the loga- 
rithmic velocity-distribution law so well as those obtained 
by Mkurai.se in his tests on sand roughness is explained by 
the fact that in our case the ratio of height of roughness 
to height of channel or pi-oc radius is considerably larger 
(approximately l/s as compared to 1/ 15 - 1/500 ). For this 
reason our profiles, particularly those in the neighborhood 
of the wall, show greater scattering since here the effect 
of the individual roughness elements is very prominent. In 
tables III and IV are given the results of the tests on the 
rough plates. The former gives for each profile the maxi- 
mum velocity U, the mean velocity of the profiles at the 
rough wall u, the kinematic viscosity v > the widths 
hi and b 2 of the smooth and rough profiles, the Reynolds 
Number U b/u , the roughness coefficient 7 #r k/v, and the 

roughness function A, and table IV, giving for each_of the 
•rough plates the nondinensional shearing stress v*/u, the 
equivalent sand roughness k Q , the value of a = k s /k, and 
the mean value of A for each plate. Of the sand-plate 
tests, only the results obtained with,, the "Hamburg sand" 
are given, since the tests with the G-bttiiigen sand are to 
be considered only as preliminary and carried out with the 
object of seeing how the results agreed with those obtained 
by JJikuradse in his tests with sand roughness. The agree- 
ment reached between his and our results is satisfactory, 
the values of A differing from those found by Nikuradse 
by 1 to 2 percent. 
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A = -& 



7 



* r 



TABLE IV 

5.75 log |, kg = equivalent sand roughness, 
a = k s /k, 5.75 log a = 8.48 - A 



Plate 
No . 



*r 
u 



cm 



= a 



XII 
III 
I 
II 
V 



VI 
IV 



XIII 
XIV 
XV 
XIX 



XX 1 1 1 
XXIV 
XXV 



XVI 
XV I 1 1 
XVII 



Sphere roughness: k=0.41 cm 



0.0689 
.0831 
.120 
.131 
.0854 



12.2 
8.92 
5 .68 
5 . 15 
9 . 65 



0. 093 
. 344 

1 . 26 
1.56 

.257 



Sphere roughness: k=0.21 cm 



0.0779 
.106 



8.9 8 
5.27 



0. 172 
.759 



Spherical segment roughness: k=0.26 cm 



0.0590 
.0631 
.076 3 
.0909 



13.8 
12.7 
9.89 
7.64 



0.031 
.049 
. 149 
. 365 



Cone roughness: k=0.375 cm 



0.065 2 
.0754 
.0894 



13.1 
10.6 
8.49 



0.059 
. 164 
. 374 



Short angle roughness: k=0.30 cm 



0.0856 
.101 
. 1 24 



8.56 
5.67 

4 . 5 3 



0.291 
.618 
1.47 



0.227 
,8 38 
3.07 
3.81 
.626 



0.819 
3.61 



0.118 
.186 
.571 

1.40 



0.159 
.437 
.996 



0.965 

2.05 

4.86 



Long angle roughness: k=0.323 cm, 0.310 cm, 0.303 cm 



XX 
XXI 
XXII 



IX 



0. 137 
.167 
.179 



4.17 

2.28 
2.33 



1.81 
3. 70 
3.56 



Sand roughness: 1:=0.135 cm 



0.0820 



7.2; 



0.222 



5.61 
11.9 
11,75 



1.64 
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The most important, from a practical consideration, 
of all the magnitudes given in our results is that of k g , 
the equivalent sand roughness, which, according to equa- 
tion (5) and figure 9, makes it possible to apply the re- 
sults to channels and pipes of other dimensions and also 
to the towed plate, (See following section.) According 
to the results obtained for the sand roughness, turbulent 
flow is fully developed for v*k s /i; > 70 and the same lim- 
it presumably holds for other types of roughness. At any 
rate our measurements all lie within the range of complete- 
ly rough or turbulent flow, since according to table III 
the value of A is practically independent of the Reynolds 
Number. In applying the results to channels and pipes of 
other heights and diameters, respectively, it is necessary, 
however, to consider the limit v*k s /\; > 70. In a pipe 

u • / T 

v* = — — and turbulent flow is obtained when 
2.83 

^/x > 198 

Some exceptional phenomena, the cause of which could 
not be explained, were observed in obtaining the velocity 
profiles for plate XII (sphere with d=4 mm, D=40 mm ) • 
To these we shall again refer in section on Velocity Dis- 
tribution at Plate XII (p. 34). Figure 16 shows a plot of 
the nondimensional velocity u/v^ r against log (y/k s ) 
obtained for all the 21 rough plates. The points all fall 
to a sufficient degree of accuracy on the straight line 

JL = 8.48 + 5.75 log J- 
v* r k s 

which was used for the determination of k g . The scatter- 
ing is very slight, showing that for all the plates inves- 
tigated the universal velocity-distribution equation is 
well satisfied. For only very small values of y/ k are 
there any systematic deviations. 

Of particular interest is the d.evendence of the re- 
sistance of rough plates, having equal roughness elements 
and similar manner of distribution, on the density^of the 
roughness. Figure 17 shows the variation of v +r /u with 
F r /F for the various plates where F r denotes the pro- 
jection of the area of all the roughness elements on a 
plane normal to tho direction of the flow and F is the 
area of the plate,, so that F r /F », 0 denotes a smooth 
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plate. With increasing F r /F the resistance at first in- 
creases as is to he expected. For the spherical rough- 
ness elements, however, the maximum resistance does not oc- 
cur at the greatest density, "but at F r /F = 0.4; that is, 
when about 40 percent of the total area is covered with 
spheres. This is also easy to understand since with the 
elements widely distributed, the entire diameter of the 
sphere is effective whereas when the packing is closer, 
only the radius or less becomes the effective absolute 
roughness. The closest packing of the spherical and spher- 
ical segment elements have a lower resistance than that of 
sand roughness whose grain size is equal to the diameter 
and height, respectively, of the spheres and spherical seg- 
ments. For both cases, according to table IV, kg/k< 1« 
The effective absolute roughness is evidently less for the 
"regular" roughness due to spheres and spherical segments 
than for the "irregular" sand roughness. 

A maximum is also found for the long angle roughness 
elements approximately at F r /F = 0.1. 



To obtain a still better understanding of the depend- 
ence of the resistance on the roughness density, we shall 
define a resistance coefficient for each type of roughness 
in the following manner. Let denote the resistance 

due to the roughness elements alone, that is, the differ- 
ence between the total resistance W of the rough plate 
and the resistance Wg of the smooth area between the 
roughness elements , 



¥ r = 77 - 



(15) 



Further, let ui c denote the velocity at a distance from 
the wall y = k equal to the height of the roughness. 
We then form a resistance coefficient 



2W 



Uf 



(16) 



p u k F 



where F r denotes as before the projected area of all the 



roughness elements on the -ple.ne normal 
the flow. (See Sadron , reference 18.) 



when the shear stress 



the wall 



to the direction of 
The velocity u k , 
p v* r 2 is known , 



may be computed from the velocity-distribution equation 
(8). With u = U for y - b 2 , we have from equation (8) 



v*r 



- 5.75 log vf- 



- 2.5 In 



y 



(1?) 
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and from the above "by integration "between the limits y - 
0 and y = b 2 , we derive the mean velocity XL of the 
velocity profile at the rough wall 

U - u = 2 • 5 v * r 

Subtracting from the above equation 



f - tt k = - 2.5 v^ T In ^~ 



we obtain 



or 



u k - u = 2.5 v, r (l * in - k -) 
* - 1 - 2 .5 ^ (in | - l) 



(18) 



For the resistance ff r , we have from equation (15) 
^r = I P V* r a - J x p v; g 2 



o r 



c r again 



U 2 v u / 1 V u / 



_!!£__ . 2 1. (1' [pa* - ?l fl!arf-l a.) 



Of - 

P u Jk F r F r 



where t' = p v 1 „ 2 denotes the shear stress at the smooth 

spaces "between the roughness elements and F x the area of 
these spaces. 

Since u/u k is known from equation (18), we have 
thus expressed the resistance coefficient^ Of in terms of 
measured magnitudes. The value of v'^g/u still depends 

somewhat on the Reynolds Number, decreasing as the former 
increases. Since the second term in the brackets. is gen- 
erally small compared to the first, we have replaced, in 
order to simplify the expression, v'wn by a mean value 
which was obtained by measurements on a smooth tinned 

v 1 

plate, the value thus found being -= — = 0.0461. The val- 
ues cf =r- and the values of Of computed therefrom by 
u 

equation (19) are given in table Y and -clotted in figure 
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18 as a function of F r /F. ?or the spherical segment's, 
cones, and angles, Of is constant for small values of 
the roughness density, so that the resistance of the rough- 
ness elements depends on the roughness density only through 
the effect of u^ . At larger values of the roughness den- 
sities the curves fall off with increasing F /F . 
The effective height of roughness "becomes smaller with the 
absolute height remaining the same. The Cf curve for 
the 4 . 1-mi llimeter spheres slopes do^n somewhat even for 
the smallest measured roughness densities; i.e., in the 
case of the spherical roughness the processes occurring a. t 
each element exert a mutual effect even at smaller rough- 
ness densities than those of the other types of roughness 
elements investigated. The t^o plates with the 2. 1-milli- 
meter spheres fall to some extent on the curve for the 
4.1-millimeter spheres, a result that is to "be expected 
from the geometric similarity of the elements. 

The three plates with the long angles behave in an 
exceptional way, the value of Cf first increasing, then 
decreasing as the roughness density increases. Tests with 
this tyoe of roughness had previously been carried out by 
Treer (reference 19) , who used values of F r /F ~ 0.5, 
0.53, 1.0, considerably larger than ours. For these large 
values of roughness densities, he finds an increase in the 
resistance with increasing roughness density, in agreement 
with our results. 

Finally, we may compare the resistance coefficients 
Cf of the roughness elements in the bounda.ry layer with 
trie usual resistance coefficients c w for unconfined flow. 
For rectangular plates of aspect ratio i/b = o° and 8/3 
(corresponding to our plates with the long and short an- 
gles), set normal to tne direction of flow c w = 2.01 and 
1.17, respectively, independent of the Reynolds Number. 
For spheres at Re = XL&fv - 2 X 10 4 which is about equal 
to the Reynolds Number of our spherical roughness, c w = 
0.47. Tnese values are plotted on figure 18 and it is 
found that they agree quite well with the Cf values at 

sma.ll roughness densities. It therefore follows that for 
small roughness densities the resistance of a roughness 
element in the friction layer is about the same as it would 
be in an unconfined flow at a velocity of flo™ eaual to 
that prevailing at the distance from the ™all y — k in 
the friction layer. For our other roughness elements we 
cannot, unfortunately, make this comparison at present 
since for these the value of c w is as yet unknown. 
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TABLE V 

■ar = roughness density, Of = resistance coefficient of 
roughness element (se< D equation (19)), u^ = velocity 
at distance y = k (see equation (18)), c Wqo = resist- 
ance coefficient in unconfirmed flow. 
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Spheres d=0.41 
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XIII 


nn q 7 

♦ UUo f 
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.480 


• Of-O 


Spher i cal segments 


XIV 


.0155 


.0039 8 




.469 


.799 


Spher i cal segments 


XV 


.0348 


.00582 




. 389 


. 767 


Spherical seginen t s 


XIX 


. 251 


.00825 




.102 


.702 


Cones 




XXIII 


.0106 


.00425 




.552 


.865 


Cones 




XXIV 


.0189 


.00569 




.561 


.832 


Cones 




XXV 


.0425 


.00799 




.463 


.790 


Short angles 




XVI 


.0151 


.00732 


1 


.20 


.75 7 


Short angles 




XVI II 


.0269 


.0102 


1 


.24 


. 691 


Short angles 




XVII 


.0605 


.0154 


1 


.19 


.607 


Long angles 




XX 


.0538 


.0188 


1.95 


.563 


Long angles 




XXI 


.0776 


.0279 


3 


.62 


.428 


Long angles 




XXII 


.152 


.0321 


2.54 


. 378 


S an d 




IX 




.00672 






.594 


Sphere , Re = 


2 x 


10 4 : c 












Rectangular plate 


, aspect ratio 


l/t> = 8/3 


• 


C Wco " 


1. 17 


Rectangular plate 


, aspect ratio 


l/"b = a> : 




C Woo = 


2.01 
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Extension of Results to Towed Plates 

The total resistance of a ship is made up, as is 
known, of skin-friction resistance, eddy-making resistance, 
and wave-making resistance - of which the first in many 
cases forms by far the greatest portion of the total. The 
frictional resistance depends very strongly on the rough- 
ness of the ship ! s surface. • The increase in the resist- 
ance due to the surface roughness may amount to 50 percent 
of the total frictional resistance. The problem of sur- 
face roughness is therefore of interest in ship design, 
particularly as regards the investigation of the resist- 
ance offered "by rough plates to tangential flow. 

Whereas, in the case of pipes and channels, the ratio 
of the absolute roughness k to the thickness of the fric- 
tion layer 8 is constant all along the pipe, the proc- 
esses occurring in the friction layer of the towed plate 
are complicated by the fact that the f r i ct i on- lay e r thick- 
ness increases from the front toward the rear so that the 
ratio k/ 8 on which the resistance depends, decreases 
from front to rear. After a short initial laminar run the 
value of " k/ 8 at the forward part of a rough plate is rel- 
atively large so fell at there we have completely turbulent 
flow. Farther hack, provided the plate is long enough and 
the roughness elements sufficiently small, there is formed 
the transitional region and "beyond this possibly a region 
of laminar flow. 

Prandtl and von Karman , the former in 1927 and 1932 
(references 13 and 14), the latter in 1921 and 1930 (ref- 
erences 12 and 4), have shown for the case of the smooth 
plate how, from the laws of flow in smooth pipes, there 
may be derived the law for smooth plates in a purely theo- 
retical way and in good agreement with the results of 
tests. Using the results of Nikuradse on pipes with sand 
roughness as a basis, the same method was applied by Prandtl 
and the author to rough plates (references 15 and 16). 
These extensions to towed plates were made possible through 
the discovery of the universal laws for turbulent flow, 
equations (7) and (8), in smooth and rough pipes. The re- 
liability of this method and its extrapolation beyond the 
range of tests conducted is assured by the universal char- 
acter of these laws. For the details of the method the 
reader is referred to the references cited. 

The diagrams obtained in this manner for the local 
resist an ce coefficients 
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C f « = 



2t_ 

PT 1 



the 



and the total resistance coefficients Of = — -g" — °^ 

pV n vl 

rough plate as a function of the Reynolds Number R = — 

and the relative roughness 1//k s (figs. 19 and 20) have 
already "been given in reference 16 cited above. It is as- 
sumed that the flow is turbulent starting from the forward 
edge. On account of the small region of laminar flow that 
is always present the total resistance is somewhat lowered, 
although for relatively long plates this lowering is prac- 
tically of no importance. 



The relations are simplest in the region of completely 
turbulent flew where the resistance coefficients depend 
only on the relative roughness. The exact law for the re- 
sistance is, however, also in this case, so c oapl i cat ed 
that it cannot he- expressed in explicit form "but only as a 
general function of a parameter (parameter z), namely, 



I = k s G(z) , W = p v 2 b k s F(z) 
(h = plate width) 



(20) 



The functions G ( z ) and F ( z ) are given in table 
VI. For practical purposes, however, it is mor^ conven- 
ient to possess simple interpolation formulas for the re- 
sistance law. For the local and total resistance coeffi- 
cients Cf 1 and Cf , we have obtained such formulas as 
approximate the accurate expression to a sufficient degree 



of accuracy and are valid for a range 
namely : 



2¥10 2 = l/k s = 10 



< 



Cf 



/ 

■ (2, 



87 + 1.58 log 



[Of .<: 



1.S9 + 1.62 log JM 

k s / 



-2 



(21a, b) 



The values of 
functions of 



Of 1 

l/k. 



smooth ulate are 



and Cf are plrtted in figure 21 as 
The corresponding formulas f<">r the 



C f • = (2 log R - 0.65 ) 



-2.3 



a 0.455 (log R) 



-2.58 



(10 6 = I 



2 10 9 ) 



(22a, b) 



Equations (21) are used for the rough plate in place ^f 
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equation (5) for the pipe. We are accordingly ahle, for 
each type of roughness whose equivalent sand roughness 
k s is known, to indicate immediately the local and total 
resistance coefficients. We shall illustrate hy means of 
an ex amp 1 e : 

Length of ship I = 150 meters, speed '*of ship v a 12 
knots per hour = 6.17 meters per second, kinematic viscos- 
ity of water at t = 15 0 G: v = 0 . 0114 centimeters per sec- 
ond, Reynolds Itfumher Re = vl / v = 8.12 x 10 8 . The total 
resistance coefficient of the smooth ship is, according to 
(22), = 1.62 x 10~ 3 . We now assume the ship through- 

out its entire length to possess a roughness of the type of 
our spherical segments (plate XIII), somewhat correspond- 
ing to rivet heads; height k = 2.6 millimeters, diameter 
d = 8 millimeters, distance apart D = 40 millimeters. 
The value found for k g (tahle V) is k s ■ 0.031 centime- 
ters. From our equation (21"b) with i/k = 4.84 x 10 5 , 
we thus compute for the rough surface C f = 2.43 x 10" 3 , 
which is a '50-percent increase as compared with the smooth 
plate. This is a very large contribution to the resist- 
ance "but the roughness assumed "by us was also considerable. 

In order to see to what extent the individual parts 
of the ship contribute tc the total resistance, we shall 
find the partial resistances W x and W 10 for the first 
and the tenth parts, respectively, cf the ten parts into 
which the ship T s length is assumed to he divided. Denot- 
ing hy W = UTi + W 2 + . . . + W 10 the resistance of the en- 
tire ship ! s surface, we find from our formula (21b ) 



! GV 

W 10 Of (I) 



I _ 1 ~ x UoV 



= 0.149 



-1° = 10 V10__./ - 0.081 

w o f ( i ) 

The first tenth portion therefore contritutes 14.9 percent 
and the last tenth only 8.1 percent to the total resist- 
ance. From this it may "be seen that the same roughness 
has a considerably more harmful effect at the hew than at 
the stern. 
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TABLE VI 

Total and local resistance coefficients Of and Of 1 and 

friction-layer thickness o/k s as functions of rela- 
tive roughness ^/^ s in turhulent flow. 

I is the length of the plate 

k g , sand roughness 

8, friction-layer thickness 

Of 1 , local resi stance coefficient 

Of, total resistance coefficient 



■ '1 

z 


I /L = 5( z ) 




F(z) 


0 f 'xlO 3 


Cf xlO 3 


fi/k q 


2 1 

3x10 ' 


1 • 615xl0 2 


1 


. 150 


9 . 84 


14.3 


l.OOxlO 1 


5xl0 2 ' 


3 . 140x 10 2 


1 


.832 


8 . 28 


11.7 


1.67X10 1 


7xl0 2 


4 . 8 32xl0 2 


2 


.490 


7 .46 


10.3 


2.33X10 1 


10 3 


7 * 5 19x 10 2 


•7 


. 45 3 


5.71 


9.18 


3. 33X10 1 


2x 10 3 


"i 7 66y10 3 

— 9 1 \J KJ /\ A. KJ 


5 


.490 


5 . 54 


7 . 35 


6.67X10 1 


SxlO 3 


2 8 7 8 x 1 0 3 


9 


.407 


4 . 99 


6 .54 


1. OOxlO 2 


5xl0 3 


5.270X10 3 


1 


.503X10 1 


4.41 


5.70 


1.67xl0 2 


7xlC 3 


7 . 841xl0 3 




. 048X10 1 


4.08 


5.22 


2.3 3xl0 2 


10 4 


1. 195xl0 4 


2 


. 843X10 1 


3. 77 


4.76 


3. 33xlO S 


2xl0 4 


2.660X10 4 


5 


. 390X10 1 


3.26 


4.05 


6.67xl0 2 


3xl0 4 


4.230X10 4 


7 


. 845X10 1 


3.01 


3.71 


l.OOxlO 3 


5xl0 4 


7 .564xl0 4 


1 


. 260xl0 2 


2.73 


3.33 


1.67xl0 3 


7xl0 4 


1.106xl0 5 


1 


. 722x10 s 


2.57 


3.11 


2.33xl0 3 


10 5 


1.654xl0 5 


2 


.400xl0 2 


2.41 


2.90 


3. 33xl0 3 


2x10 s 


3.581x10 s 


4 


.581xl0 2 


2.15 


2.56 


6.67X10 3 


3x10 s 


3.611xl0 5 


! 6 


. 693xl0 2 


2.01 


2.39 


l.OOxlO 4 


5xl0 5 


9 .888x10 s 


1 


.080x10 s 


1.86 


2.18 


1.67xl0 4 
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Formulas (21) are valid only provided the. absolute 
roughness may be considered the same all along the ship T s 
length. The practical problem sometimes arises of finding 
the resistance of a plate whose roughness varies at differ- 
ent regions. For this case, too, our formulas are still 
applicable though in passing from one roughness to the 
other, it is necessary to modify the procedure in the spe- 
cial manner described "briefly "below. 

We shall consider a plate that is of roughness fc g 
over a portion of length l 1 and roughness kg for the 
remaining distance 1 2 • The resistance of the first por- 
tion may be found as usual from formula. (21b). In comput- 
ing the resistance of the remaining portion, however, it 
is incorrect to .take x = I x as the initial coordinate 

but another coordinate x = 1 1 1 different from I x must 
be taken. Since in the transition region "between the two 
roughnesses the momentum at the "boundary layer varies con- 
tinuously, l x l is determined, from the condition 

"°2 ° 1 

i.e., lx f is that fictitious distance from the leading 
edge at which, for the roughness k S2 the loss in momen- 
tum in the friction layer would be the same as that for 
roughness k s at distance l 1 . To determine l x ? we 
must use equation (20) from which 

h- = &(*!>, hi = Hn*t 

The first equation gives the parameter z 1 at the end of 
the first portion of the plate. The relation "between z x 
and z x T follows from (23), namely, 



W 



i = p v 2 | k Gi F(z x ) = p v 2 d k s (f 2 X ») 



k 



no that . F(zi') = r~ F(2i) (24) 

It is thus possible to compute z 1 J from z 1 and there- 
fore also to compute 1 l T . The parameter z 2 at the end 
cf the second portion of the plate is then obtained from 

&( Z2 

K s 2 
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The resistance of the second portion of the plate is 
W 2 - o v 2 h k Sg [F(z 2 ) - f (l x O] 

so that the * re si stance of the entire plate is 
ffx + W 2 - p v 2 h k S2 F(z a ) 

We shall again illustrate "by an example: Let the ship 150 
meters long, considered in the previous example, have its 
roughness for the first 50 meters increased to three times 
as great as "before, k s = 0.093 centimeter; so that 

|l ■ 50 m, k g = 0.093 cm, l a = 100 m, fe^ = 0.031 cm 

Frcm li /k g = G(z x ) = 5.38 x 10 4 there is obtained "by 
(211)): 

fif (% x ) = 3.53 x 10" 3 

From the tahle for G(z) and F ( z ) , we find hy interpola- 
tion 

* 3.7? x 10 4 : F(zJ = 95.5 
and from equation (24), we have : 
F( z 1 I ) = 286.5 , z i ' — 1. 22xl0 5 

&(z i) = -1- = 2.06x10 , l x l - 63.9 m 
*s 2 

— « 5.29x10 = G(z 2 ): z 2 = 2.36X10 5 



k 



S 2 



k Sa F(z 2 ) 

F(z 2 ) = 641, C f (li+la) = 2 l ^ — = 2.64x10 3 

The increase in the resistance due to the greater rough- 
ness at the forward third portion of the ship's surface 
thus amounts to 9 percent. 

The computation that was here carried out for a plate 
with two different roughnesses may of course he easily 
extended to three or more roughnesses. 

We thus possess all the necessary data for computing 
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the required resistance for any rough surface of a ship 
whose 11 equivalent sand roughness" k g has "been obtained 
by experiment in our channel. 

Our test apparatus, which may now "be considered as 
having attained its final form, makes it possible to in- 
vestigate every type of roughness that may occur in prac- 
tice, particularly the surface roughness of ships. It is 
only necessary by the above process to determine in each 

/ v *k\ 

case the roughness function A = A [ )« In the turbu- 

\ ** / v*- 
lent region, where A is constant and independent of — -r* 

it is possible, with the aid of the equivalent sand rough- 
ness k s , according to formulas (5) and. (21a, b), to extend 
the computations to other pipes and channels as well as to 
tewed plates. In this way the plate towing terts for the 
determination of the ship's resistance, which a- e often 
very difficult to carry out, are replaced "by the much sim- 
pler tests in our channel. 

For the region of transition where the roughness 

y . V- 

function A defends on the extension of the results 

v 

to other pipes and channels and to towed plates is prssible 

/ v *k\ 

only if the roughness function A = A I— ^- ) is accurately 

known. This function probably depends also on the quality 
of the roughness, and further investigation is necpssary. 

Velocity Distribution at Plate XII 

In measuring the velocity field in the neighborhood 
of the rough wall of plate XII (diameter of sphere = 4 mm, 
separating distance = 40 mm (see fig. 6a)), a peculiar 
phenomenon was observed which appears to he so contrary to 
all cur present-day views on turbulent velocity distribu- 
tion that it seems worth while to describe it here "briefly* 

In figures 22 and 23 are shown the velocity distribu- 
tion diagrams for two planes normal to the rough plate, 
the lines connecting points of eaual velocity. Figure 22 
is for the plane immediately behind the last row of spheres 
while figure 23 is for a plane 40 millimeters behind this 
row. In the former the sphere a lying directly in front 
of the pitot tube naturally causes large velocity losses^ 
From an examination of both cf the figures, however, it is 
clear that the smallest velocities occur at the free spaces 
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"between the spheres and the greatest velocities occur "be- 
hind the rov/s of spheres :;here a greater reduction in ve- 
locity should "be expected. In figure 23, for example, the 
spheres b and c, 4 centimeters from the measuring plane, 
and sphere e, 8 centimeters away, clearly show excess 
velocities and the same is true for spheres To and c in" 
figure 22. The shielding effect of the rows of spheres is 
therefore, so tc speak, negative. This phenomenon was ob- 
served for this one plate alone having the maximum sphere 
separation. The phenomenon clearly shows that the disturb 
ance produced by the individual spheres is so strong that 
it reaches' over to the neighboring spheres, so that the 
processes at each of the spheres have an effect on each 
other even at the relatively great distance of separation 
with D/d = 10. The effect also expressed itself in the 
value of the resistance coefficients C f for the individ- 
ual roughness elements (fig. 18). 

In order to check whether this phenomenon was not due 
to the limited cross-sectional area of the channel flow 
and possibly caused by the interference effect of the op- 
positely lying wall, the same roughness was enlarged so as 
to be geometrically similar (sphere diameter = 10 mm, dis- 
tance apart = 100 mm) and the p^ate, 150 centimeters long, 
was investigated in the large Gottingen wind tunnel (jet 
diameter = 2.25 meters). Exactly the same phenomenon was 
again observed, showing that it was not brought about by 
the fact that the stream was restricted to a finite cross 
section but that it was a pure fricrt ion-layer phenomenon. 
No satisfactory explanation of this phenomenon has as yet 
been found. Possibly complicated secondary flows have 
something to do with it. At any rate, it appears that the 
turbiilent mixing processes behind obstacles in the neigh- 
borhood of walls are in some respects very different 
from those behind obstacles in an unconfined flow. 



SUMMARY 



Based on the universal laws of turbulent velocity dis 
tribution at rough and smooth walls, there is in the pres- 
ent work presented a method that allows surface roughness 
tests and in particular, measurements on the roughness of 
ship surfaces" to be carried out in a much simpler manner 
than was done heretofore. The types of roughness investi- 
gated were in the form of flat, rough plates installed in 
a square-section rectangular channel, the other three wall 
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always being smooth. Twenty- one plates of various rough- 
ness were Investigated, the roughness elements "being the 
following: spheres of diameter 0.41 and 0.21, respective- 
ly, spherical segments, cones, and "short" and "long" an- 
gles. The ah so lute roughness or roughness height k was 
2 to 4 centimeters. The pattern, according to which the . 
elements were distributed, was the same for each plate 
while the density of the roughness was varied. 

For each rough plate the nonsymmetrical velocity dis- 
tribution was measured for six different velocities at the 
outlet section of the channel and the pressure drop ob- 
tained by orifices in the rough wall. Preliminary tests 
established the fact that the velocity profiles at the smooth 
and rough trail J exerted no interference effect on each oth- 
er, each velocity distribution being the same as if the 
walls were all smooth or all rough, respectively. From the 
measured nonsymmetrical velocity distribution, it is there- 
fore possible with the aid of the v el o ci t y- di s t r ibut i on 
law, to determine the sheer stress Tg at the smooth wall 
and the required shear stress at the rough wall T r is 
then found from the equation 

b d*o 

1 & p dx 

where b is the channel height. For almost all the rough- 
nesses and velocities investigated, the resistance was in- 
dependent of the Reynolds Number. In order to apply the 
results of the measurements conveniently to channels and 
pipes of other hydraulic radii T-u and to towed plates, 
there was determined fer each of the roughnesses investi- 
gated an 11 equivalent sand roughness" fc g which is the 
grain size of the sand roughness as used in the tests of 
Nikuradse, and has the same resistance as the correspond- 
ing roughness elements. 

Besides depending on the relative roughness r^/k, 
the resistance also depends on a second roughness parame- 
ter, namely, the roughness density F r /F where F r is 
the total projected area of the roughness elements on a 
plane normal to the direction of flow and F is the plate 
area. In determining the dependence of the resistance on 
the roughness density F r /F the fact was established that 
the maximum resistance does not occur at the maximum rough- 
ness density but at a considerably lower value. Further- 
more, for each of the rough plates was determined the co- 
efficient Of of the roughness clement referred to the 
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velocity u^ at distance y = k and defined "by the equa- 
tion 

°f = f-~" 

P u k F r 

where W r denoted the resistance of the rough portion of 
the surface (after subtracting that of the smooth areas 
"between the roughness elements). It was found that for 
nearly all the types of roughness elements the value of 
Cf is independent of the roughness density F r /F for 
small values of the former and falls off sharply for large 
values. The resistance coefficient Cf of the individual 
roughness elements practically agrees in value with the. 
coefficient in an unconfined stream for the same Reynolds 
Number . 

For each of the elements whose equivalent roughness 
k s was determined in the manner indicated, the results 
could k extended immediately to find the resistance of 
the towed plate with the aid of a diagram given "by L. 
Prandtl and H. Schlichting in a previous work on the re- 
sistance of rough plates. For the region in which the tur- 
bulent flow-resistance law applies, simple interpolation 
formulas are given for the local and total resistance co- 
efficients of the towed plate. 



Translation "by S. 
Hat ional Advi s ory 
for Aeronautics. 



Reiss, 
Commi 1 1 ee 
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TABLE III 

U= maximum velocity, u = mean velocity of the profile at the 

1 dp 

rough wall, - — == pressure drop, V = kinematic viscosity 
p dx 

p ss density, = distance of maximum velocity from smooth 

wall, "b a = distance of maximum velocity from rough wall, 

v o. = / — ^ = "shear stress velocity 11 at smooth wall, v - = 
*g J p * r 

— - shear stress velocity at rough wall, t = shear 
P 

stress at smooth wall, T r = shear stress at rough wall, 

u-L v A -k 

Re = — = Reynolds Number, r = roughness coefficient. 
v v 
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O.OI 16 


5 


658 


523 


1 165 


0,0Il6 


6 


778 


626 


1622 


0,0Il8 


i 


313 


240 


260 


0,01 14 


2 


384 


297 


403 


0,0II5 


3 


500 


384 


696 


0,0Il6 


4 


568 


438 


908 


0,01 17 


5 


646 


499 


1184 


0,01 13 


6 


746 


572 


1585 


0,01 I 2 


i 


3ii 


259 


186 


0,01 16 


2 


385 


322 


284 


0,OII5 


3 


498 


417 


473 


0,0114 


4 


585 ! 


489 


648 


0,01 13 


5 


662 


554 


833 


0,01 14 


6 


809 


678 


1247 


0,0II3 








Continued 



1,51 

1,49 
1,49 
J, 52 
J. 52 

1,06 
1,16 
1.27 
1,26 

1.42 
1,26 

0,91 
0,96 
1,01 
1 ,01 
0,96 
0,96 

0,98 
0,90 
0,90 
0,88 
0,90 
0,92 

1,16 
I,l6 
I,l6 
1,16 
1,16 
I,l6 



2,48 
2,50 
2.50 
2-47 
2-47 

2,93 
2,83 
2,72 
2.73 
2,57 
2,73 

3,05 
3,oo 

2.95 
2,95 
3,oo 
3,00 

2,90 
2,98 
2,98 
3,oo 
2,98 
2,96 

2,52 
2,52 
2,52 
2,52 
2,52 
2,52 



= 0,41 
19,4 
23,3 
28,8 

33,i 
39,4 

23,5 
29.1 
37.2 
42,2 

52.3 
60,4 

29,6 
36.6 

48,5 
54>o 
62,8 
74,2 

3i,4 

38,7 
50,2 

57.° 
65,0 

75,o 
22,2 

27»5 
35.6 
41,8 
47.3 
57,8 



14.1 


0,073 


109 


6,80 


16,2 


0,069 


132 


8,24 


1 8,8 


0,069 


164 


10,2 


21,7 


J 0,068 


190 


II. 8 


27,0 


0,067 


224 


14,0 


I5,i 


0,086 


i 108 


8,34 


17,1 


0,086 


122 


10,5 


22,4 


0,087 


175 


J3,4 


24.5 


0,088 


202 


15,5' 


29,9 


0,090 


249 


19,0 


34-3 


0,085 


285 


21,7 


13. 1 


0.124 


107 


10,6 


15.3 


0,122 


131 


12,95 


22,6 


0,120 


168 


16,6 


25,4 


0,118 


218 


19,1 

• 


28,0 


0,1 19 


225 


22,2 


33,i 


0,1 16 


261 


25.9 


i3,7 


0,130 


I0 7 


ii,3 


*5>7 


0,129 


I30 


13,8 


20,2 


0,132 


167 


17,8 


23,7 


0,131 


187 


20,0 


26,3 


0,128 


222 


23,6 


29,2 


0,128 


258 


27,5 


14,8 


0,087 


86 1 


7,85 


17,7 


0,085 


123 


9,81 


22,9 


0,084 


161 


12,8 


26,2 


0,086 


190 


15,2 


29,8 


0,085 1 


214 j 


17,0 


35,5 


0,085 1 


263 1 


21,0 


pages. 









1 2, 1 6 
12,21 
12,19 
12,22 
12,14 

8,91 
8,95 
8,94 
8,93 
8,93 
8,91 

5,68 

5,69 
5,66 
5,68 
5,68 
5,7o 

5.20 
5,18 
5,13 
5,M 
5,M 
5,ii 

9,6i 
9,69 
9,64 
9.64 
9,67 
9,63 
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mate 

No. 


Profil< 
No. 




u 


idp 

odx 


V 




h 






u 


7 ! 

0 

»o 1 i 


— * . io -2 

V 


A 






cm/sec | cm/sec 


cm/sec 2 


cm 2 /sec 


cm 


cm 


cm/sec 






£>l 1 












Sphere 


roughness: 


d = k — 0,21 


cm 










VI 


i 


316 


269 


i54 


0,01 16 


i,36 


2,63 


20,8 


13,9 


0,077 


I09 


3>76 


8,96 




2 


39o 


333 


234 


0,01 16 


1,36 


2,63 


25.7 


17,6 


0,075 


134 


4,65 


9,02 




3 


491 


419 


365 


0,01 16 


i.35 


2,64 


3M 


20,2 


0,077 


169 


5,86 


8,99 




4 


566 


481 


480 


0,01 14 


1,35 


2,64 


37»3 


22,8 


0,078 


198 


6,87 


8.97 




5 


664 


56/ 


656 


0,01 14 


r,35 


2,64 


43,8 


27,7 


0,078 


232 


8,07 


8,99 






806 


660 


966 


0,0118 


i.35 


2,64 


53.^ 


33.2 


' 0,082 


2 73 


9,46 


8.95 


IV 


I 


3^5 


261 


242 


0,0120 


1,06 


2,91 


27,7 


13.6 


0,108 


108 


AM 


5,24 




2 


396 


319 


34 8 


0,0121 


1,07 


2,90 


33,8 


18,1 


0,105 


I3O 


5,87 


5.31 




3 


572 


459 


704 


0,0122 


1,1 1 


2,86 


48,8 


22,6 


i 0,107 


186 


8,40 


5.32 




4 


646 


5M 


888 


0,0125 


I, II 


2,86 


55,i 


25,9 


0,107 


205 


9,26 


5,25 




5 


75i 


606 


1 196 


0,0121 


1,01 


2,96 


64,1 


30,1 


0,105 


246 


1 1 , 1 


5.24 



Spherical segment roughness! A =0,26 cm 



I 


3H 


279 


124 


0,0110 


1,61 


2,38 


16,4 


14.5 


0,061 


II 4 


3,89 


13.80 


2 


389 


345 


177 


0,0110 


1,66 


2,33 


20,3 


I7>8 


0,059 


I 4 I 


4,82 


13,85 


3 


495 


438 


278 


0,0112 


1,76 


2,23 


25,9 


20,7 


0,060 


I76 


6,02 


13,86 


4 


574 


508 


367 


0,0112 


1,62 


2,37 


30,0 


23.5 


0,060 


220 


6,98 


13.78 


5 


U4O 


573 


462 


0,0112 


1,63 


2,36 


33,8 


28,2 


0,058 


23I 


7,00 


J 3»77 


5 


°J U 


/jj 


750 


0,0110 


1,79 


2,20 


43.3 


37.1 




T 

3 UA 


TO 1 




1 


380 


334 


181 


0,0113 


1,56 


2,43 


20,9 


17,8 


0,063 


134 


4,83 


I2J6 


2 


497 


43i 


298 


0,0114 


1,62 


2,37 


27.3 


21,6 


0,065 


174 


6,26 


12,62 


3 


572 


498 


39o 


0,0116 


1,62 


2,37 


3L5 


23,8 


0,064 


197 


7,08 


12,67 


4 


700 


609 


570 


0,0114 


1,57 


2,42 


38,5 


32,0 


0,062 


246 


8,82 


12,70 


5 


834 


73o 


800 


0,0114 


1,62 


2,37 


45.9 


38,4 


0,059 


292 


IO,5 


12,75 


1 


382 


328 


223 


0,0113 


1,54 


2,44 


25,0 


16,5 


0,078 


135 


5,77 


9,85 


2 


5°2 


428 


375 


0,0118 


1,52 


2,46 


32,8 


20,4 


0,078 


169 


7,58 


9,77 


3 


504 


483 


472 


0,0118 


1,58 


2,40 


36,9 


23,8 


0,075 


190 


8,16 


9,86 


4 


687 


592 


695 


0,0117 


1,62 


2,36 


45.o 


30,3 


0,073 


244 


10,0 


9,94 


5 


817 


710 


977 


0,0114 


1,61 


2,37 


53.4 


34,8 


0,074 


284 


12,3 


10,04 


1 


3i6 


264 


198 


0,0114 




2,68 


23,8 


i4>i 


0,091 


105 


5,40 


7,68 


2 


386 


319 


294 


0,0113 


1,27 


2,58 


29,2 


17,4 


0,091 


132 


6,74 


7,66 


3 


480 


401 


452 


0,0116 




2,68 


36,2 


20,4 


0,090 


159 


8,16 


7,66 


4 


563 


47i 


620 


0,0117 


M7 


2,68 


42,5 


26,6 


0,088 


184 


9,49 


7,66 


J 


671 


553 


880 


0,0115 


1,25 


2,60 


50,6 


28,5 


0,093 


224 




7,62 


6 


818 


673 


1305 


0,0112 


1,27 


2,58 


61,8 


34,5 


0,092 


281 


14.4 


7.58 








Cone roughness! k = 


0,375 cm 










1 


321 


281 


' 130 


0,0110 


1,56 


2,43 


18,2 


^4,4 


0,064 


116 


6,21 


13.07 


2 


386 


335 


185 


0,0109 


1,70 


2,29 


21,9 


16,7 


0,066 


141 


7.54 


13.08 


3 


488 


427 


291 


0,0109 


1,64 


2,35 


27.7 


21,2 


0,064 


178 


9,54 


13,10 


4 


574 


497 


4°3 


0,0108 


1,72 


2,27 


32,6 


24,6 


0,067 


212 


n,3 


13,04 


5 


668 


578 


542 


0,0106 


i,74 


2,25 


38,0 


27,6 


0,067 


252 


13,4 


13.08 


1 


307 


260 


147 


0,0108 


i,49 


2,49 


20,1 


13.4 


0,08l 


113 


6,99 


10,52 


2 


384 


326 


218 


0,0108 


i,49 


2,49 


25,2 


16,4 


0,08l 


143 


8,74 


10,50 


3 


495 


423 


350 


0,0107 


1,48 


2,5° 


32,5 


20,5 


0,077 


184 


«/4 


10,59 


4 


567 


485 


452 


0,0108 


i,49 


2,49 


37,2 


24,0 


0,075 


207 


12,9 


10,59 


5 


662 


564 


605 


0,0106 


1,52 


2,46 


43.4 


28,0 


0,074 


248 


15.4 


10,56 


6 


890 


756 


1076 


0,0109 


1,50 


2,48 


58,4 


36,2 


0,076 


325 


20,1 


io,54 
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i late 

No. 1 


.riOUK 

No. 


- U 
cm/sec 


w 

cm/sec 


idp 
gdx 
cm/sec 2 


V 

cm 2 /sec 


h 

cm 


L 

°2 

cm 


cm/sec 


cm/sec 


V 9 r 

U 


b 

br I 


IO 2 

V 


A 


XXV 


i 


3io 


257 


176 


0,0107 


1 

i,33 


2,62 i 


23,4 


12,2 i 


0,096 


115 


8,19 


8,41 




2 


388 


324 


267 1 


0,0107 


1,32 


2,63 


29,2 


16,9 


0,090 


144 


10,3 


8,48 




3 


476 


397 


401 


0,0106 


1,36 


2,59 


35,9 


20,7 | 


0,089 1 


172 


12,7 


8,51 




4 


564 


470 


562 


0,0107 


i,35 


2,60 


42,5 


23,3 


0,091 


208 


M,9 


8,50 




5 


668 


555 


788 


0,0106 


1,42 


2,53 


50,4 


277 


0,089 


249 j 


17,8 


8,54 




6 


787 


660 


1086 


0,0105 


1,30 


2,65 


59,3 ' 


32,0 


0,089 


296 


21,2 


8,54 






" Short angle" 


roughness I k = 


0,30 cm 










XVI 


1 


3iS 


267 


175 


0,0112 


1,52 


2,48 


23,1 


12,5 


0,089 


113 


6,24 


8,54 




2 


389 


328 


254 


0,0112 


1,46 


2,54 


28,2 


15,5 


0,088 


139 


7,62 


8,52 




3 


497 


422 


415 


0,0110 


1,40 


2,60 


36,i 


19,3 


0,086 


181 


9,9i 


8,54 




4 


566 


482 


534 


0,0110 


1,40 


2,60 


4M 


23,0 


0,083 


206 


n,3 


8,58 




5 


666 


563 


737 


0,0111 


1,46 


2,54 


48,4 


27>5 


0,085 


240 


13,2 


£ C-7 




6 


789 


671 


1030 


0,0108 


1,46 


2,54 


57.2 


32 Q 


O.O84 


259 


16,0 


8,64 


XVIII 


i 


380 


310 


300 


0,01 14 


1,24 


2,76 


3i,6 


14,6 


0,102 


133 


8,38 


6,6> 




2 


4Q8 


411 


5H 


0,0117 


1,21 


2,79 


4M 


18,4 


0,100 


170 


10,7 


6,66 




3 


S62 


463 


654 


0,0116 


1,21 


2,79 


46,8 


20,0 


0,102 


194 


12,2 


6,66 




4 


666 


547 


921 


0,01 16 


1,21 


2,79 


55,4 


23,8 


0,102 


230 


M,4 


6,66 




5 


796 


661 


1310 


0,0114 


1,22 


2,78 


66,2 


31.3 


O, IOO 


279 


i7,5 


6,76 


XVII 


i 


313 


239 


301 


0,0111 


1,07 


2,92 


3i,6 


13,1 


0,138 


113 


8,60 


4,39 




2 


l8s 


296 


440 


0,0113 


1,07 


2,92 


38,9 


l6,0 


0,138 


136 


10,4 


4.43 




3 


*\OI 


400 


722 


0,01 16 


1,1 1 


2,88 


5°«6 


18,4 


0,123 


172 


'3,2 


4>63 




4 


567 


45i 


925 


0,0116 


1,1 1 


2,88 


57,3 


21,6 


0,124 


195 


14,9 


4,^2 




5 


674 


527 


1292 


0,01 11 


1,1 I 


2,88 


68,1 


28,4 


0,127 


242 


18,5 


4,53 




6 


773 


616 


1693 


0,0113 


1,01 


2,98 


78,1 


31,2 


0,124 


272 


20,9 


4-59 






"Long angle M 


roughnessl k - 


0,32 cm 




9, 16 




XX 


i 


3i6 


240 


334 


0,01 16 


0,72 


3,i8 


32,9 


10,8 


0,134 


10b 


4,20 




2 


388 


291 


500 


0,01 17 


0,80 


3,io 


40,4 


12,2 


0,140 


126 


11,1 


4> I( > 




3 


515 


399 


865 


0,0120 


0,72 


3,18 


53,6 


13,9 


0,133 


t At 
107 


r 4'4 


4» J 9 




4 


57 8 


445 


1092 


0,0120 


0,74 


3>i6 


60,2 


20,1 


0,135 


l88 


16,2 


4 iS 




5 


Ai- Q 
650 


498 


1408 


0,01 17 


0,75 


3,i5 


68,5 


21,6 


0,142 


2 20 


18,9 


j. 1 > 
-t> 1 - 
















k = 


:o,3i cm 








XXI 


I 


313 


214 


450 


0,0117 


0,66 


3,30 


38,8 


ii,3 


0,192 


I 06 


10,3 


2,2 2 




2 


388 


270 


648 


0,0117 


0,62 


3,34 


48,1 


13.3 


0,l82 


131 


I2 ,7 


2,25 




3 


488 


348 


975 


0,0117 


0,58 


3,38 


60,6 


15,6 


0,1 72 


T A- 


tA t 
1 0, 1 


2,29 




4 


566 


413 


1281 


0,0116 


0,63 


3,33 


7<>>3 


17,9 


0,167 


193 


I 0,0 


2,3° 




5 


722 


529 


2022 


o.oi 14 


0,68 


3,28 


89,7 


28,4 


0,162 


2 ^O 


24,4 


2,36 
















k = 


0,30 cm 








XXII 


i 


ore 


214 


392 


0,0109 


0,69 


3,27 


38,7 


9.8 


0,182 


1 *5 


IO,8 


2.30 




2 


^Q4 


270 


605 


0,0108 


0,67 


3,29 


48,4 


12,3 


0,180 


145 


13.6 


2.33 




3 




339 


920 


0,01 12 


0,59 


3.37 


59,o 


14.4 


0,171 


171 


l6,0 


2,41 




4 


^67 


389 


1270 


0,OII2 


0,65 


3.31 


69,6 


16,1 


O.I79 


200 


18,8 


2,32 




5 


666 


456 


1740 


0,0108 


0,64 


3,32 


81,8 


19,7 


0,180 


244 


22,9 


2,30 




6 


749 


512 


2200 


O.OI IO 


0,64 


'3,32 


92,0 


24,0 


0,180 


270 


25,3 


2.31 








Roughne 3 s f Hamburg sands } < 


= 0,135 cm 








IX 


i 


319 


270 


181 


0,0116 


1,25 


' 2,62 


22,0 


14.5 


0,083 


106 


2,56 


7,27 




2 


420 


354 


297 


0,0118 


1,22 


2,65 


29,0 


I9,i 


0,081 


138 


3,32 


7>'9 




3 


501 


423 


416 


0,0120 


1,19 


2,68 


34,6 


21,9 


0,080 


162 


3,89 


7,20 




4 


564 


472 


516 


0,01 16 


1,28 


2,59 


38,9 


23,3 


0,083 


188 


4,53 


7,46 




5 


654 


546 


692 


0,0119 


1,18 


2,69 


45,i 


28,1 


0,082 


212 


5»i2 


7»« 




6 


763 


639 


933 


0,0121 


1,16 


2,71 


52,6 


31,3 


0,082 


244 


5,86 


7,ii 




log Re — — 

Figure 1.- The resistance factor X=-2_ £ of rough pipe* as a function of the Reynolds ^ 

^ P5 £ 

Number Re=^ and the relative roughness r/k for sand roughness (after Nikuradse) . ^ 



Figure 2.- Test set-up for the roughness 
measurements. 

wk=r water tank; zr= inlet pipe; ak= chan- 
nel for initial run; mk= measuring chan- 
nel; gm= velocity measuring apparatus; 
kp=* centrifugal pump; vk= storage channel; 
mb=^ measuring tank; fl = pulley; s bivalve 
between pun.p and tank; a bivalve between 
tank and inlet pipe. 





Figure 3.- Teat set-up for surface 
roughness measurements. 
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Figs. 4,5,8 



■3200- 



Figure 4.- Test chan- 
nel for 
surface roughness 
measurement. 
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Figure 5.- Dimensions and dis- 
tribution of roughness elements 
I, II = spheres, 
k=0*41 and 0.21 cm.; 
Ill* spherical segments, k= 0.26 
cm.; IV = cones, k -0.375 cm; 
V=* short angles, k=*0.30 cm; 
VI = long angles, k=0.30 cm; 
0.31; 0.32 cm. 



o il= W0 cm/sec • 'cm/sec • u- 525 cm/sec • u= 550 cm/sec 
• -¥50 " • '500 n 6 -550 " 

//////////////// //////////////////// ////////////////////////////////< 



'ft* 




-770mm- 



Figure 8.- Lines of equal velocity at outlet section of channel 
obtained with hoth walls smooth. 
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Fig. 6a— i 





Figure 6a- Spheres .Plate XII. Figure 6b- Spheres, Plate III, 




Fig. 6c- Spheres, Plate I. Fig.6d- Spheres f Plate II. Fig.6e- Spheres f Plate V, 



d = 0.41 cm: k - 0,41 cm; 
D — 1,0 cm; *, = 126 cm. 



d = 0.41 cm; k = 0,41 cm; 
D = 0,6 cm; *j = 1.56 cm. 



d — 0.41 cm; k = 0.41 cm: 
dichteste Packung: k $ = 0.257 cm. 





W f .f * **• * * » * * ** * < 
* , * * i • * f * » » » ♦ ' 



Figure 6f- Spheres .Plate VI 4 




Figure 6g- Spheres .Plate IT. 

d = 0,21 cm; k = 0.21 cm: 
D = 0,5 cm; k. = 0.759 cm - 




Figure 6h- Spherical segments, 
Plate XIII. 



it = o.S cm; * =» 0.26 cm: 
D = 4.0 cm; k s = 0.031 cm. 



Figure 61- Spherical segments, 
Plate XI?. 

d = o,8 cm: k = 0.26 cm; 
D = 3.0 cm; k 5 = 0,049 cm. 



Figure 6(a~i).- Roughness elements of the plates investigated. d*diam. of 

roughness, D=mean distance of separation of elements; 
k=» height of roughness; k 8 = equivalent sand roughness. 
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Fig.6k--s 




Fig. 6m,. Cones .Plate XXIII. Fig.6nr Cone^l^^Sl€ Fig. 60.- Cones .Plate XX?. 




Figure 6(k~s).- Roughness element* of the plates investigated. d«diam. of 

roughness. 15 -mean distance of separation of elements; 
k = height of roughness; k g - equivalent sand roughness. 
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Figs. 7, 9 
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Figure 7.- Nonsymmetrical velocity distribution in channel 

having smooth and rough walls. Roughness is that 
of Hamburg sand, k * 0.135 cm. 
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Figure 9.- The resistance factor A for completely rough flow with sand 
roughness as a function of the relative roughness %/k 0 . 
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Fig. 10 
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figure 10.- The nondiaensional velocity distribution for the 
spherical roughness, d = k= 0.41 em* 
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Pig. U 




N.A.C.A. Technical Memorandum Ho. 823 



Fig. 12 
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Figure 15.- The nondimeneional velocity distribution for the long 
angle reughneee, k =0.32. =0.31; = 0.30 cm. 
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rig. 14 




0,5 0,6 

Figure 14.- The nondimensional velocity distribution for the 
short angle roughness, k = 0.30 cm* 
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Figs. 17,18 
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Fi-rure 17.- The dependence of the 

resistance on the 
roughness density fJfV \ 



— — as a function of — for 
pu 2 * 
sphere , spherical segment , cone , 
and an^-le roughness. 
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Figure 18.- The resistance coefficient for the 
indivdual rougmess element 

2W r F r 
Cf - : — as a function of the roughness density 
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Figure 19.- The local resistance coefficient C 1 ** 



for rou^rh sand 



plates as a function of the Reynolds Number 2- and 

V 



the relative roughness 



, (k = sand roughness) 
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Figure 20.- The total resistance coefficient ?f = — g — for rough sand plates as a 

pv - 

function of the Reynolds Number and the relative roughness 

(k s = sand roughness) . 
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Fig. 21 
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Figure 21.- The total and local resistance coefficients 

and C- for completely rough flow as a 
function of the relative roughness : 
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Figs ,22,23 
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Figure 23.- Lines of equal velocity in plane at distance 10 d=4.0 cm 
behind last row of spheres (Plate XII). 



